FAST 工程 主动 反射 面 健康 监测 系统 中 温度 信号 多 元 回归 估 
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摘 要 : 大 型 结构 健康 监测 系统 中 ,温度 测 点 故障 会 造成 安全 隐患 。 基 于 FAST 工程 主动 


反射 面 健康 监测 系统 ,提取 并 分 析 了 9 个 测 点 数据 间 的 线性 相关 性 , 将 候选 变量 G 
最 优 回 归 子 集 ， 建 立 多 元 线性 回归 模型 ， 融 合 正常 测 点 数据 对 故障 测 点 进行 估 测 。 针 对 变 

间 多 重 共 线 性 问题 进一步 应 用 了 岭 回归 法 ， 选 定 岭 参数 为 6。 使 用 F 检验 与 拟 合 度 检验 
型 显著 性 和 有 效 性 , 使 用 不 同时 日 数据 验证 了 模型 的 佑 测 精 度 。 结 果 证 明 ， 多 元 线性 回归 模 
型 比 单 变量 具有 更 高 的 拟 合 度 与 精度 ， 均 方 根 误差 为 0.475C， 岭 回归 法 稳定 性 更 高 ， 均 方 
根 误差 仅 为 0. 435°C。 
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结构 健康 监测 是 通过 布设 大 量 传感器 , 监测 并 评估 结构 当前 状态 的 一 种 技术 , 在 工程 结 
构 的 安 安全 监测 中 有 着 广泛 应 用 。500 米 口径 球面 射电 望远镜 FAST(Five- hundred-meter 
Aperture Spherical Telescope) 大 型 射电 望远镜 为 口径 500 米 具有 主动 反射 面 的 大 型 钢 结构 工 
程 。 文 献 [2] 讨 论 了 FAST 工 程 结构 受 力 复杂 性 ， 由 于 其 结构 特殊 、 空 间 跨 度 大 且 环 节 众多 ， 
在 观测 过 程 中 主动 反射 面 受 控 变 形 致 结构 受 力 复杂 , 为 保证 其 安全 工作 , 建设 了 主动 反射 下 
健康 监测 系统 ， 通 过 监测 关键 结构 的 应 力 与 环境 信息 ， 实 现 FAST 望 远 镜 的 健康 状态 评估 。 
文献 [3] 指 出 , 在 索 网 施工 过 程 中 , 该 系统 所 监测 圈 梁 实时 应 力 最 大 约 60Mpa， 始终 小 于 设计 
安全 值 201. 5Mpa， 保 持 在 安全 范围 内 ， 并 指出 了 温度 信号 对 应 力 分 析 的 重要 性 。 
在 FAST 运 行 期 间 ， 工 作 过 程 中 索 网 主动 变形 会 使 结构 受 力 更 为 复杂 ， 以 2019 年 7 月 份 监 
测 数据 为 例 ， 所 有 轿 梁 格 构 柱 测 点 应 力 值 最 大 值 为 137. 37Mpa， 最 小 值 为 -135. 88Mpa， 仍 处 
于 安全 范围 内 ， 但 部 分 测 点 应 力 变 化 范围 较 大 ， 应 力 变 化 最 大 的 测 点 其 极 差 为 167. 17Mpa。 
在 望远镜 后 续 工 作 中 , 应 力 的 监测 对 结构 实时 健康 状态 评估 、 长 期 疲劳 损伤 评估 等 仍 具有 重 
要 意义 。 
大 型 工程 结构 空间 尺寸 大 , 不 同 部 位 温 升 不 同 , 温差 与 温度 效应 会 使 结构 应 力 发 生 较 大 
变化 ， 文 献 [ 和 指出 FAST 工 程 结构 中 ， 温 度 荷载 对 圈 深 结构 的 刚度 起 控制 作用 ， 分 布 测 点 温 
度 信息 是 评估 结构 状态 的 重要 依据 ”"。 且 在 FAST 工 程 中 ， 应 用 了 416 只 光纤 Bragg 光 栅 (Fiber 
Bragg Grating, FBG) 应 变 传感器 "监测 应 力 ， 共 包含 圈 梁 及 格 构 柱 应 力 测 点 100 只 ， 主 索索 力 
测 点 316 只 ， 虽 测 点 分 布 位 置 不 同 ， 但 所 有 测 点 均 需 根据 不 同 测 点 的 实际 温度 补偿 应 变 传 感 
器 的 温度 应 变 交叉 敏感 问题 ， 文 献 [7] 进 一 步 讨 论 了 分 离 结构 应 力 的 温度 效应 对 状态 评估 的 
必要 性 ， 温 度 信息 缺失 将 直接 导致 测 点 应 变 信息 不 可 靠 。 
实际 工程 现场 环境 复杂 ， 大 量 传感器 长 期 工作 , 传感器 及 数据 通道 存在 一 定 故 障 率 。 检 
修 发 现 , FAST 主 动 反射 面 健康 监测 系统 已 有 个 别 光 纤 光 栅 解 调 仪 出 现 故 障 , 且 存在 数 只 传 感 
器 数据 异常 ， 已 确定 异常 测 点 数目 占 比 约 为 8%， 使 部 分 监测 数据 无 法 正常 获取 。 但 受 限 于 现 
场 施 工 条 件 ， 以 及 部 分 测 点 安装 位 置 特殊 ， 故 障 传感器 常 无 法 第 一 时 间 修 复 ， 数 据 的 缺失 将 
导致 健康 监测 系统 性 能 下 降 ， 形 成 安全 隐患 。 
使 用 临近 测 点 值 代替 故障 测 点 温度 信息 ， 就 FPBG 传 感 器 而 言 ， 因 其 常 串 联 布设 , 若 数 据 
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通道 发 生 损坏 , 会 使 临近 的 一 组 传感器 均 失 效 。 使 用 较 远 传感器 监测 值 代替 , 存在 距离 过 远 ， 


光照 影响 温差 较 大 ,代替 不 准确 等 问题 , 这 也 是 每 测 点 都 布设 温度 传感器 的 原因 。 实 现 故障 
测 点 的 温度 信号 准确 估 测 ， 对 提高 健康 监测 系统 可 靠 性 有 着 重要 意义 。 
1 传感器 布设 与 数据 提取 分 析 

FAST 主 动 反 射 面 结构 健康 监测 系统 圈 梁 格 构 柱 关 键 应 变 测 点 的 温度 信息 ， 使 用 FBG 温 


度 传 感 器 监测 。 传 感 器 工作 原理 是 外 界 温 度 引 起 光栅 周期 以 及 有 效 折射 率 发 生变 化 ,使 反射 
波长 偏 移 ， 通 过 测量 中 心 波长 变化 量 ， 获 取 传 感 器 所 测 温 度 。 
为 有 效 探究 远 距离 布设 的 温度 测 点 之 间 的 相关 性 ， 提 取 FAST 边 缘 圈 梁 支 承 格 构 柱 对 应 


测 点 的 温度 信息 ， 每 个 被 监测 的 格 构 柱 提取 一 路 测 点 信息 ，FAST 工 程 圈 梁 格 构 柱 共 50 个 ， 
有 10 个 格 构 柱 装 有 测 点 , 被 测 格 构 柱 编号 与 测 点 分 布 如 图 1 所 示 , 传感器 安装 于 如 图 2 所 示 轿 
梁 支 座 处 的 靠 内 侧 的 水 平 拉 杆 中 的 绿色 杆 件 ， 测 点 之 间 最 远 距 离 为 500m。 
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图 1 测 点 分 布 示意 区 图 2 传感器 安装 位 置 
Fig.1 Distribution map of measuring points Fig.2 Sensor mounting position 
监测 系统 数据 采集 频率 为 IHz， 为 降低 数据 量 ， 滤 除 高 频 干 扰 ， 且 避免 由 于 数据 长 度 不 


统一 造成 的 分 析 困 难 。 对 原始 数据 进行 平均 处 理 ， 平 均 时 窗 为 10min， 取 每 10 分 钟 内 温度 数 
据 平 均值 作为 该 时 刻 温度 数据 ， 每 测 点 每 日 得 144 个 数据 样本 。 提 取 2019 年 7 月 份 监测 数据 ， 


对 应 记 作 Ty ~ Tig, ， 其 中 31# 格 构 柱 测 点 因 采 集 通道 检修 关闭 无 数据 记录 ， 共 提取 9 个 测 点 
六 息 。 因 环境 温度 循环 周期 往往 为 县 夜 ， 按 天 对 数据 进行 分 组 ，7 月 16 日 数据 如 图 3 所 示 。 
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图 3 温度 监测 数据 
Fig.3 Temperature monitoring data 
可 见 不 同 测 点 之 间 温 差 较 大 ， 最 大 温差 超过 5'C， 波 动 规律 也 有 一 定 区 别 ， 但 其 趋势 大 
体 一 致 。 两 数据 样本 分 别 记 作 关 -[x.2,,..x,], Y 2Lly. y» y,lo VSEASISEEEAS H1 ZR 
性 相关 性 ， 相 关系 数 公 式 为 : 


Y -30,-» 
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bem C (1) 
(5e YOU 
RP: n 为 样本 长 度 ，x%,y 为 各 自 均 值 ， 得 不 同 测 点 温度 数据 相关 系数 矩阵 如 表 1， 可 见 


不 同 测 点 之 间 存 在 高 度 线性 相关 ， 相 关系 数 最 低 为 0.959。 
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Tab.1 Correlation coefficient matrix of temperature data 

Ti Te- Tug o Cam Ta Ta fe Ea Ta 
Tu 1000 0.973 0.986 0.969 0.982 0.984 0.997 0.996 0.996 
Te 0973 1.000 0.980 0.981 0.974 0.984 0.967 0.969 0.982 
Tı 0.986 0.980 1.000 0.991 0.993 0.988 0.981 0.975 0.985 
Tis 0.969 0.981 0.991 1.000 0.994 0.992 0.959 0.959 0.972 
Taw 0.982 0.974 0.993 0.994 1.000 0.995 0.974 0.973 0.981 
Tə 0.984 0.984 0.988 0.992 0.995 1.000 0.977 0.980 0.988 
Tae 0.997 0.967 0.981 0.959 0.974 0.977 1.000 0.996 0.995 
Tas 0.996 0.969 0.975 0.959 0.973 0.980 0.996 1.000 0.997 
Tas 0.996 0.982 0.985 0.972 0.981 0.988 0.995 0.997 1.000 


2 线性 回归 模型 


2.1 模型 简介 


当 变 量 之 间 高 度 线性 相关 时 ， 常 用 线性 回归 方法 建 模 ， 自 变量 个 数 为 1 时 ， 记 为 x*， 因 
量 记 作 y， 其 一 般 模 型 记 为 : 


H 
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ar 


in| 


y-f,tBxce (2) 
RP: B» ARERIO € 为 随机 误差 。 
健康 监测 温度 估 测 应 用 中 , 受 光 照 不 均匀 的 影响 , 一 元 回归 在 某 些 时 刻 预 测 精度 较 差 ”。 
引入 更 多 变量 ， 纤 合 不 同 测 点 的 温度 影响 ， 可 取得 更 优 的 预测 结果 。 当 变量 个 数 为 多 个 
Xp AD Xp, 称 为 多 元 线性 回归 ， 模型 记 为 : 


y-fyt Boat Bx, +p, X +E (3) 
RP: Bp Dee L NRG pLdetkAj3S SIRENA, PRA E EEE X: 
y=XB+E (4) 
yı Po E 
E 
Y= Y» B= P, [NL (5) 
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AF: YA, XEFE, BABET AI. EWIE, KEWL 
误差 项 符合 正 态 分 布 6; ~ N(0,a ),i=1,2,…,n 。 售 去 误差 项 的 影响 ， 最 小 二 乘法 求解 模型 
的 自 变量 系数 ， 系 数 向 量 估计 为 : 


B-(X'X)"X'Y (7) 
利用 所 得 自 变量 估计 系数 值 B, , ,., 记 , 即 可 建立 经 验 回归 方程 : 


y- fy* B, x... B,x, (8) 


式 中 : x m 加 为 自 变 量 输入 值 ，y 为 模型 输出 估 测 值 。 
对 提取 的 9 个 测 点 数据 ， 取 其 一 测 点 作为 因 变 量 即 被 估 测 点 ， 其 余 作为 自 变 量 输入 ， 即 


可 通过 训练 建立 多 元 线性 回归 模型 ， 当 被 估 测 点 数据 缺失 时 ， 可 利用 自 变量 输入 ,对 其 信号 
进行 估 测 。 


22 最 优 自 变量 选取 
将 1# 格 构 柱 测 点 作为 被 估 测 点 ， 其 余 8 个 测 点 温度 信息 作为 候选 变量 ， 其 模型 变量 的 不 
同 分 组 子 集 有 28 -1 种 。 利 用 16 日 数据 训练 模型 ， 将 变量 按照 相关 性 降序 排列 逐个 增加 ， 划 
分 为 8 组 。 
结合 C, 统计 量 和 R2 自由 度 调整 复 决定 系数 两 选取 准则 来 确定 最 优 变 量 的 子 集 : 


RSS, 
C,-(n-m Jo R n+2p (9) 


m 


AP: n=144 为 样本 长 度 ，m=8 为 全 部 待 选 变 量 个 数 ，p 为 模型 所 选 变 量 个 数 ，RSS,, 为 选择 
全 部 变量 建 模 的 残 差 平 方 和 和 ，RSS, 为 模型 残 差 平方 和 。 


R;-1- "l alg (10) 
n—-p-1 
式 中 : R 为 复 决 定 系数 即 模 型 拟 合 度 
pjese (11) 
TSS 
式 中 : 75S5 为 因 变 量 y=T, 离 差 平方 和 。 
RSS = D0 —»,) (2 
i-l 
TSS = (y, - »y (13) 
i-l 


RP: AREE, y 为 模型 预测 值 ，y 为 实际 值 均 值 。 


变量 分 组 与 对 应 模型 C, A Rz 数值 如 表 2 所 示 ， 选 取 C, fU HR; 最 大 ， 即 为 最 优 变量 
子 集 ， 可 见 八 元 作为 变量 时 模型 拟 合 度 最 优 。 
表 2 变量 分 组 与 模型 C, 、R; 数值 
Tab. 2 Variable grouping and C,, R; value of model 


Groups Variables C p R 

Univariate Tsor 130.51 0.9937 
Two variables  T36#, Tai 98.4908 — 0.9944 
Three variables  T36#,T41#,T46+ 56.04 0.9954 
Four variables Tis Tos Tai Taos 20.95 0.9963 
Five variables — Tii Tos Tac Tato, Tas 22.0 | 0.9963 
Six variables Tus Toi Tos T36, Tato, Tas 22.0 | 0.9963 
Seven variables — Tg Tii Tous T26#, Tos, Tat, T4 13.29 0.9965 


3 ”多重 共 线性 问题 与 岭 回归 模型 


由 表 1 知 ， 所 有 测 点 两 两 之 间 均 存在 极 强 相 关 性 ， 对 于 多 元 线性 回归 模型 ， 会 因 多 重 共 
线性 使 模型 稳定 性 出 现 问题 ， 即 对 模型 可 能 存在 一 组 数 操 ,后 …, 使 : 
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此 


时 设计 和 矩阵 秩 ramk(X)< 书 +1， 意 味 着 公式 (7) 中 | XTX [e 0 ， 而 模型 系数 向 量 估 
EME p AEE D) =P (XTX ， 此 时 方差 阵 对 角 元 素 过 大 ， 使 自 变量 系数 向 量 B 的 
估计 精度 下 降 ， 模 型 对 输入 数据 的 


ko +kix; +. +k, x m0, -L2..n) ( 


pi 
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14) 


小 变化 会 过 于 敏感 ， 预 测 结果 不 确 


定性 增加 ， 模 型 的 稳 


定性 和 精确 度 降 低 。 
岭 回 归 法 针对 多 元 线性 回归 模型 中 的 多 重 共 线性 问题 进行 了 优化 ,向 量 系数 岭 估 计 定 
义 为 : 


XTX + 人 17 满 秩 可 首 ， 接 近 奇 异 程度 比 原 矩阵 减 小 。X "7 X + AI 随 着 | 
型 方差 减 小 ,但 向 量 系数 估计 值 p 将 偏离 原 值 ， 使 模型 偏差 增 


元 线性 


4 


BA) - OC X AD? XTY (15) 
AP: 4 >0 称 为 岭 参 数 。 可 见 给 X7X SUNT —4 1E UA 


E 阵 L2 范 数 惩罚 项 ， 保 证 


可 归 模型 ， 在 应 用 时 ， 需 合理 选取 岭 参 数 。 


系数 趋 近 于 0， 说 明 岭 回归 法 训练 期 间 ， 对 自 变 量 有 筛选 作用 。 
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Estimation of vector coefficients 
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Fig.4 Ridge parameter map 


4 ”模型 检验 与 评价 
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元 线性 回归 模型 建立 后 ， 需 判断 模型 是 否 可 靠 ， 通 过 F 检 验 确定 模型 建立 的 回归 方程 


叭 参数 增 大 而 增 大 ， 模 
大 ，4=0 模 型 退化 为 普通 多 


HIRIRA, 将 8 个 自 变量 回归 系数 对 应 岭 迹 绘制 如 图 4 所 示 , 按 原 则 选取 使 所 有 向 量 系 
的 最 小 的 一 个 岭 参 数值 ， 确 定 4 = 6 ， 建 立 岭 回归 模型 。 由 图 
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显著 性 ， 统 计量 公式 为 : 


式 中 : RSS 为 模型 残 差 平 方向，ESS 为 回归 平方 和 和 ，n=144 为 档 


ESS/p 
RSS /(n— p —1) 


(16) 


ESS = Y,» a7) 
i-l 


vb: y 为 模型 预测 值 ，y 为 实际 值 均值 。 


模型 拟 合 度 尺 2 检验 ， 见 公式 (11) ， 值 越 接近 于 1 说 明 


以 Ti# 数 据 作为 被 估 测 点 ， 其 他 测 点 数据 作为 输入 变量 。 使 用 16 


并 对 17 


^ 18H 


数据 进行 估 测 检验 ， 使 用 均 方 根 误差 RMSE 评 价 模型 


说 明 模 型 预测 效果 越 好 : 


4.2 测试 结果 


RMSE-|LY0,-»? 8) 
ni 


模型 拟 合 度 越 高 。 
日 数据 训练 建立 模型 ， 


的 估 测 效果 ， 


本 长 度 ，p=8 为 所 选 变量 个 数 。 


值 越 小 
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对 照 前 述 分 组 ,建立 从 一 元 到 八 元 不 同类 型 的 线性 回归 与 岭 回 归 模 型 ， 并 对 模型 进行 估 
测 效 果 评 价 ， 检 验 及 模型 估 测 结果 见 表 3。F 检 验 统计 量 最 低 为 4972.29， 依 旧 远 大 于 查 表 所 
得 Fa(1,135)= 3.911， 所 有 模型 通过 显著 性 检验 。 拟 合 度 均 在 0.99 以 上 ， 满 足 拟 合 度 要 求 。 随 
着 融合 变量 的 增加 ， 模 型 测试 的 RMSE 也 随 之 减 小 。 
表 3 模型 检验 与 应 用 测试 结果 
Tab. 3 Model test and application results 


16th model fitting 17th estimate — 18th estimate Estimated mean value 
Groups F test R Fitting degree | RMSE/'C RMSE/'C RMSE/'C RMSE/'C 
Univariate 22517.92 0.9937 0.2 0.68 0.66 0.67 
Two variables 12784.23 0.9945 0.19 0.71 0.56 0.635 
Three variables 10411.94 0.9955 0.17 0.52 0.46 0.49 
Four variables 9593.81 0.9964 0.15 0.52 0.44 0.48 
Five variables 7665.63 0.9964 0.15 0.52 0.45 0.485 
Six variables 6878.33 0.9967 0.15 0.53 0.43 0.48 
Seven variables 5876.86 0.9967 0.15 0.53 0.43 0.48 
Eight variables 5418.98 0.9969 0.14 0.52 0.43 0.475 
Ridge regression 4972.29 0.9966 0.1 0.46 0.41 0.435 
其 中 ， 一 AE. RAAR RA, ORDRURTEE RC 
y 2 HEN 
ELTER, ASTRITIN S, i 信 回 归 法 模型 精度 更 高 。 
==- Univariate 
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图 5 真实 值 与 模型 预测 值 对 比 
Fig.5 Comparison of true value and model predicted value 


5 结 语 


提取 FAST 工 程 主动 反射 面 健康 监测 系统 中 跨度 500m 范 围 布 设 的 9 个 测 点 温度 监测 数 
dm. 分 析 知 两 两 间 存 在 极 强 的 线性 相关 性 ,将 1# 测 点 作为 被 估 测 点 , 研究 了 融合 其 余 测 点 数 
据 对 其 估 测 的 方法 ， 并 对 比 了 不 同 变 量 集 时 模型 的 估 测 效果 。 

多 元 线性 回归 对 比 一 元 的 回归 模型 , 无 论 是 在 拟 合 度 还 是 实际 估 测 应 用 中 , 都 具有 更 好 
的 效果 ， 当 输入 变量 增加 至 三 元 时 ， 测 试 中 RMSE<0.5C 。 而 岭 回归 对 比 原始 多 元 线性 回归 ， 
模型 拟 合 度 略 有 降低 ， 但 可 有 效 避 免 变量 间 的 多 重 共 线 性 问题 ， 对 比 17、18 日 数据 的 测试 ， 
其 预测 精度 要 更 优 于 多 元 回归 ， 算 法 具有 更 强 的 抗 干扰 能 力 ， 稳定 性 好 ， 两 日 平均 RMSE 仅 
730.435'C 。 

结果 表明 , 当 某 测 点 出 现 故 障 时 ,可 使 用 岭 回归 法 利用 近期 历史 数据 训练 并 建立 多 元 线 
性 回归 模型 ， 将 正常 测 点 温度 信息 作为 输入 ， 对 故障 测 点 输出 值 进行 估 测 。 方 法 具有 较 高 精 
EE, 可 用 于 补充 健康 监测 系统 缺失 的 节点 温度 信息 , 在 故障 测 点 维修 期 间 保持 健康 监测 系统 
的 正常 工作 。 但 该 方法 使 用 近期 日 期 数据 进行 拟 合 预测 ， EATA AAIE 
案 ， 若 测 点 故障 过 久 ， 可 用 数据 时 日 相差 较 大 ， 由 于 日 照 天 气 差别 等 原因 ， 设 备 温 度 场 变化 
规律 迁移 ， 估 测 精度 将 较 低 ， 长 期 数据 缺失 的 佑 测 方法 仍 需 进一步 研究 。 
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Abstract: The temperature measurement point malfunction in the large-scale structural 
health monitoring system could cause safety hazards. The paper extracts 9 temperature 
measurement points’ data of the FAST project active reflective surface health monitoring system 
and analyzes the linear correlation of them. In order to establish a multiple linear regression model, 
the variables are grouped and the optimal subset is selected. The fault measuring point is estimated 
by fusing the normal measuring points' data. Focusing to the multiple collinear issue between 
variables, the ridge regression method is applied and the ridge parameter is set to 6. Using F-test 
and fit-degree, the significance and validity of models are tested. Based on different time data, the 
estimation accuracy is verified. The results show the multiple linear regression model has a higher 
degree of fit and accuracy than one variable model, considering its RMSE is 0.475 C. Besides, the 
ridge regression model is more stable, considering its RMSE is 0.435 C. 

Keywords: Five-hundred-meter Aperture Spherical Telescope(FAST); temperature sensor; 


health monitoring system; multiple linear regression; ridge regression 


